Nucleation threshold and deactivation mechanisms of nanoscopic cavitation nuclei 
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The acoustic nucleation threshold for bubbles trapped in cavities has theoretically been predicted 
within the crevice theory by Atchley & Prosperetti [J. Acoust. Soc. Am. 86, 1065-1084 (1989)]. 
Here, we determine this threshold experimentally, by applying a single pressure pulse to bubbles 
trapped in cylindrical nanoscopic pits ("artificial crevices") with radii down to 50 nm. By decreasing 
the minimum pressure stepwise, we observe the threshold for which the bubbles start to nucleate. 
The experimental results are quantitatively in excellent agreement with the theoretical predictions 
of Atchley & Prosperetti. In addition, we provide the mechanism which explains the deactivation 
of cavitation nuclei: gas diffusion together with an aspherical bubble collapse. Finally, we present 
superhydrophobic nuclei which cannot be deactivated, unless with a high-speed liquid jet directed 
into the pit. 



I. INTRODUCTION 



Water can be ruptured at much smaller tensile stresses 
than theoretically expected [1]. The reason for this dis- 
crepancy is the existence of small inhomogeneities in the 
liquid, which exist even when special care on the clean- 
liness of the water has been taken [2|]. The inhomo- 
geneities, whatever their origin might be, have received 
the generic name "cavitation nuclei", while the bubble 
generation produced in this way is termed heterogeneous 
nucleation. Cavitation nuclei are generally long-lived and 
it is believed that they consist at least in part of a vol- 
ume of gas [3]. This observation excludes the possibil- 
ity of the inhomogeneities being free spherical gas bub- 
bles, as these are unstable. To account for stable gaseous 
cavitation nuclei, two types are distinguished in the lit- 
erature: bubbles stabilized by a skin (see Ref. Q and 
references therein) and bubbles trapped inside a surface 
defect (" crevice model" ) 0, i, H, 0, i, H • 

The principle of the crevice model dates back to 
1944 and has found extensive qualitative experimen- 
tal evidence over the years. Greenspan and Tschiegg [10], 
for example, reported that removing particles larger than 
0.2/im in diameter increased the tensile strength of wa- 
ter to about 200 bar (see also Refs. [2, 11]). Others 
found that the addition of suspended particles lowers the 
nucleation threshold [H, [H, [11, [H, [11, [13, while pre- 
experimental pressurization of water increases the nucle- 
ation threshold[5",T8']. Although these findings are in line 
with the general idea of the crevice model, none of the 
experiments could quantitatively verify the theoretical 
crevice model as developed by Atchley & Prosperetti Q 
in 1989. One of the reasons is that their predictions are 
valid for a single cavity of a well-defined shape, while in 
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practice the liquid usually contains a wide variety of nu- 
clei of different sizes and shapes. Even in ultrapure water 
with a controlled number of microparticles, the sizes of 
the nuclei present on the microparticles can exhibit size 
variations yielding a wide range of thresholds p^ . 

A step forward was achieved by Bremond et al. who 
were able to create monodisperse cavitation nuclei by 
trapping gas inside cylindrical holes of well-defined shape 
etched in silicon surfaces using standard lithography 
techniques [l^, [2O, [2l|. Not only was the position and 
size of the nuclei perfectly controlled, but the nucleation 
event itself was also highly reproducible, so that it could 
be followed in time with stroboscopic methods without 
the need of expensive high-speed cameras. 

Both conditions, the reproducibility of the experiment 
and the monodispersity of the nuclei present at fixed po- 
sitions, are important ingredients of this paper. Here, 
we have downscaled the micropits of Bremond et al. two 
orders of magnitude, so that it becomes possible to ex- 
perimentally verify the theoretical predictions made in 
the framework of the crevice model. This is the first aim 
of this paper. Secondly, we explore the mechanisms lead- 
ing to the deactivation of nuclei after a single nucleation 
event. In addition, we show that superhydrophobic cav- 
itation nuclei can nucleate hundreds of times, without 
being deactivated. Our observations and interpretations 
have implications for an increased understanding of the 
behavior of cavitation nuclei down to length scales of a 
few tens of nm ("surface nanobubbles" ) [22]. 



II. BRIEF THEORETICAL DESCRIPTION 

A comprehensive description and development of the 
crevice model can be found in the paper of Atchley & 
Prosperetti [8] with extensions to any axisymmetric ge- 
ometry, including cylindrical cavities, given by Chappell 
& Payne In both papers, the authors predict the 

various nucleation thresholds as function of crevice shape, 
gas tension, and receding contact angle. Atchley & Pros- 
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peretti argue that a true nucleation event must be the 
result of the loss of mechanical stability of the nucleus, 
and calculated the threshold for two situations: the first 
threshold denotes the pressure at which the nucleus in- 
side a crevice begins its unstable growth, and is the low- 
est pressure value belonging to a bubble reaching either 
the critical radius of curvature Rc^ or the receding radius 
of curvature Rr^ i.e. the radius of curvature at which the 
receding contact angle Or is reached. The second thresh- 
old is equivalent in definition, but holds for the bubble 
growing outside of the crevice mouth. The lower value 
of the first and second nucleation threshold is the one 
for which the bubble grows explosively out of the cavity. 
For the case of a bubble trapped in a cavity with volume 
Vc with its interface at the crevice mouth, the (second) 
nucleation threshold is given by 

^"-^ R -^""^ Vc^{7r/3)g{0)R^' 

with PL the liquid pressure, py the vapor pressure, Vq 
the initial volume of the gas, Pg^o the initial gas pres- 
sure in the bubble, a the liquid-gas surface tension, and 
{7T/3)g{0)R^ the volume of the spherical-cap-shaped bub- 
ble with radius R as it expands above the cavity [8] (see 
Fig. HD. Here, g{0) = (2 + (2 + sin^ l9) cos (9) is a geo- 
metric function depending on the contact angle 6. The 
right-hand side of Eq. [1] represents the expanding forces 
caused by the vapor and gas pressure, py -\- Pg^ respec- 
tively, while the left-hand side represents the collapsing 
forces due to the liquid pressure and surface tension, re- 
spectively. Eq. [1] implies that, for nucleation to occur, 
the expanding forces should exceed the collapsing forces 
(condition 1). Secondly, this condition should persist for 
increasing i?, i.e. d{pg -\-py)/dR > d{pL + 2a/R)/dR 
(condition 2). 

Let us now consider the case - which we will examine 
experimentally - of a cylindrical cavity with radius Tc and 
depth dc, under the assumption that the initial gas-liquid 
interface at the crevice mouth is approximately flat (i.e. 
we assume a negligible effect of the hydrostatic pressure 
and gas-saturated water), so that we can write Vq = Vc = 
Tvr'^dc. For Or > 7r/2 the nucleation threshold is the 
pressure needed to pull the bubble beyond its minimum 
radius while it expands from the cavity, i.e. i^min = Rr = 
Tel s'mOR. Now, Eq. [T]can be rewritten as 
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FIG. 1: Cylindrical cavity with its dimensions. The initial 
gas- liquid interface is flat (dashed line), while the expanding 
bubble has a radius of curvature R and contact angle with 
the flat surface. 
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FIG. 2: Graph of the expanding forces (right-hand side of 
Eq. [T] dashed line) and the collapsing forces (left-hand side 
of Eq. [T] solid lines) for a bubble expanding from a cylindri- 
cal pit with rc — 250 nm, dc — 500 nm, pg^o — 10^ Pa — Pv, 
Pv = 73-10^ Pa and Or — 100°. The expanding forces (dashed 
line) show two branches corresponding to the possible solu- 
tions of R. The upper branch reflects the initially flat bubble 
during its first expansion phase: R decreases from R — Roo to 
R = Rr; the lower branch shows the solutions for the bubble 
expanding during its second phase: R increases from R = Rr 
to larger sizes. The collapsing forces are shown for two cases: 
1) PL = -0.20 MPa (thin solid line) and 2) pL = -0.486 MPa 
(thick solid line). For case 1, the expanding forces dominate 
over the collapsing forces and the bubble will expand, until at 
R = 0.49/im a stable equilibrium is reached. No nucleation 
will occur in this case. In case 2 the expanding forces are 
larger than the collapsing forces for all possible solutions of 
R, and as a result the bubble will grow explosively. 



^ 3p^,o 2(7 sin Or 

" ^ 3 + {r^dMOn)/ sin^ Or n ' 

This prediction will be verified experimentally in this 
paper. At the point R = Rmin = Rr the bubble experi- 
ences the maximum collapsing force. Any further reduc- 
tion of the liquid pressure will make the bubble expand, 
thus reducing the surface tension pressure, and the bub- 
ble will grow explosively with the contact angle fixed at 
Or. Equations [1] and [2] are illustrated in Fig. O for an 
air bubble trapped in a cylindrical pit with dimensions 



Tc = 250 nm and dc = 500 nm and withp^^o = 10^ Pa— p^, 
Py = 73- 10^ Pa and Or = 100°. The graph shows the 
expanding and the collapsing forces for two cases: 1) 
PL = -0.20 MPa (Blake threshold for a free bub- 

ble with Ro = Vc) and 2) pL = -0.486 MPa (predic- 
tion of Eq. [2]). From this plot it is readily seen that the 
expanding pressure in the second case is always larger 
than the collapsing pressure, and that d{pg -\-py)/dR > 
d{pL^2a/R)/dR. 

For Or < 7r/2 the nucleation threshold is much more 



complicated to calculate. At R = Rmin = the col- 
lapsing force due to surface tension is indeed maximum, 
but now d{pg -\- py)/dR < d{pL + 2a/R)/dR and there- 
fore PLiRmin) cauuot be the nucleation threshold. Also 
Pl{Rr) is not the correct threshold, as this gives a sta- 
ble equilibrium in the upper branch of the compressive 
force curve. Instead, the threshold needs to be found by 
numerically solving d{pg -\- py)/dR = d{pL + 2a / R) / dR^ 
which is an implicit equation as the contact angle de- 
pends on the radius R, through sin 6 = Vc/R^ as long as 
the receding contact angle has not yet been reached. 

In the prediction oi pl the gas term is significant for 
cavities down to a few hundred nm in radius and will 
therefore be taken into account in the present analysis. 
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FIG. 3: Sketch of the experimental setup. 



III. MATERIALS 8^ METHODS 
A. Experimental setup 

The experimental setup to investigate the nucleation 
behavior of bubbles trapped in well-defined cavities is 
similar to that used by Bremond et al. [l^, [iO, Uli 
and sketched in Fig. [3l Cavitation is induced by a fo- 
cused shock wave generator (Piezoson 100, Richard Wolf 
GmbH) consisting of piezoelectric elements mounted on 
a spherical cap at the bottom of the liquid bath, which is 
filled with 1 liter of water (Milli-Q Synthesis AlO, Milli- 
pore). The cavitation activity is recorded optically with 
a CCD camera (Flowmaster, LaVision) through a long- 
distance microscope (Model K2, Infinity). Illumination 
is provided by a fiash lamp in refiection mode. The liquid 
pressure PL is obtained with the help of a calibrated glass 
fiber hydrophone (FOPH 500, RP Acoustics). The pres- 
sure is derived by measuring the refiected intensity of the 
laser beam at the fiber tip, which depends on the density 
of the water as affected by the local pressure [26]. At 
the acoustic focus the pressure signal is typically charac- 
terized by a pressure peak (duration ~ 1 jis) followed 
by a negative pressure phase (~ bfis). The intensity 
of the pressure pulse can be varied in twenty discrete 
steps. Since the smallest possible pressure decrease at 
the acoustic focus min(pi,) = — 3.2MPa is already too 
large for our purpose, the samples are translated hori- 
zontally (away from the acoustic focus) along the line of 
sight, until the pressure signal is sufficiently weak that 
nucleation does not occur at the smallest pressure drop, 
but only at larger pressure decreases. The correspond- 
ing pressure signals are recorded using a low pass filter 
and averaged over 25 recordings to reduce the noise. A 
typical recording of the pressure signal obtained 25 mm 
out of focus is shown in Fig. 2] and corresponds to the 
experiment with the sample containing pits of 246 nm in 
radius. 
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FIG. 4: Three pressure signals with increasing strength 
recorded at the chip surface 25 mm out of focus, correspond- 
ing to the experiment with sample B {tc — 246 nm). Each 
line is the mean of 25 recordings. From these signals the min- 
imum pressure can be extracted: -0.24 MPa (thin solid line), 
-0.35 MPa (dashed line) and -0.54 MPa (thick solid line) 



B. Samples with nanopits 

The substrates of interest are silicon pieces of 5 x 5 mm^ 
diced from a Si(lOO) wafer. The nanoscopic cylindrical 
pits are directly etched into the substrate by a focused ion 
beam (FIB) in a 6 x 6 square pattern, with 200 jiia dis- 
tance between the pits. The resulting 1x1 mm^ pattern 
is located at the center of the chip. In order to facili- 
tate detailed imaging by atomic force microscopy (AFM) 
and scanning electron microscopy (SEM) identical pits 
are etched the chip corner. Four samples (A-D) were 
studied with the following dimensions (determined with 
SEM) of the nanopits: A) Vc = 495 nm; B) Tc = 246 nm; 
C) Tc = 53 nm; D) Tc = 50 — 60 nm. In sample A-C 
the pattern consisted of uniformly sized pits, with depth 
dc = 2rc, while in sample D each column of pits had 
different depths (75, 100, 200, 300, 500, and 1000 nm), 
which inffuenced the radial pit size per column by a few 
nm (50, 50, 50, 55, 57, and 60 nm resp.); see Fig.[5]for the 
corresponding SEM pictures. After production the sam- 
ples were cleaned ultrasonically in ethanol (15 minutes). 
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FIG. 6: Cavitation bubbles nucleated from cylindrical pits 
with radius Tc = 495 nm (left), Tc = 246 nm (middle) and 
rc = 53 nm (right), and depth dc = 2rc for a pressure pulse 
with pm -3.2 MPa. 



FIG. 5: SEM images of samples A-D, respectively. The bars 
in each picture denote 500 nm. The ellipsoidal shapes for 
sample A and B are due to drift. 



followed by an oxygen plasma (5 min.), a chemical clean- 
ing step using a fresh (5:1) Piranha mixture (30 min.) 
and again an ultrasonic bath in ethanol (15 min.). This 
yielded clean and completely wetting substrates, which 
were characterized by a smoothly dewetting contact line 
(if not, the whole process was repeated). Subsequently, 
the samples were hydrophobized with 1-H ^2H- 
perfluorodecyltrichlorosilane following Ref. . The ad- 
vancing and receding contact angles on the surface were 
Oa = 124° and Or = 100°. After immersion in water it 
was confirmed with AFM in tapping mode that a hori- 
zontal gas-liquid meniscus was present at the mouth of 
the pits. 



IV. RESULTS 
A. Nucleation from gas-filled nanopits 

To see whether it is possible to nucleate bubbles from 
nanopits as small as 53 nm in radius, samples A-C were 
immersed in the liquid bath allowing air to be entrapped 
in the pits. In successive experiments the samples were 
put at the acoustic focus of the shock wave generator 
and subjected to a pressure pulse with pm = m.m{pL) = 
—3.2 MPa. This value is sufficiently below the nucle- 
ation thresholds of the three samples, i.e. -0.23 MPa, 
-0.48 MPa, and -2.59 MPa respectively (see Eq.[2]), to ex- 
pect nucleation of bubbles from the nanopits. The cam- 
era and flash were triggered a few /is after passage of the 
negative part of the shock wave to capture the expand- 
ing bubbles at maximum sizes. The result is depicted in 
Fig.[6]for the samples A-C, respectively. Sample A and B 



showed almost perfect bubble patterns, with each bubble 
corresponding to the position of the nanoscopic cavita- 
tion nucleus. In each experiment, the cavitation nuclei 
had to be 're-activated' (filled with air again), since it was 
not possible to nucleate bubbles a second time without 
taking the sample out of the water first [20] . With sample 
C a maximum amount of 34 bubbles could be nucleated in 
the first experiment, implying that it is indeed possible to 
nucleate bubbles from such small cavities. While sample 
A and B showed perfect reproducibility, the number of 
bubbles nucleating from sample C declined dramatically 
in later experiments, even when the negative pressure 
amplitude was increased to -7 MPa. Presumably, small 
contaminant molecules had decreased the contact angle 
locally, leading to completely wetted nanopits. To test 
this possibility, the old hydrophobic coating was stripped 
off with an oxygen plasma and the sample was hydropho- 
bized again through the cleaning and coating steps de- 
scribed before. This process indeed re- activated part of 
the nuclei (~ 80% of the pits) though the number of 
bubbles declined again in successive experiments. 

B. Determination of the experimental nucleation 
threshold 

The experimental nucleation threshold of the nano- 
sized cavities can be obtained by moving the samples to a 
position in the liquid bath where the pressure drop Pm is 
sufficiently small that no nucleation occurs. By lowering 
Pm stepwise the cavities will nucleate at a certain nega- 
tive pressure amplitude which is the experimental nucle- 
ation threshold. To observe the bubbles optically, it is 
not sufficient for the negative pressure to be low enough, 
but it should also last long enough in time. A lower limit 
to the time At the bubble needs to grow to visible size 
is estimated, by first estimating the critical size i?c,o to 
be optically observable. We take i?c,o ^ 3 image pixels 
= 3 pix X 2.9 /im /pix = 8.7 /im. Now, using 21L \2m 

/ \l/2 

^ = ( i ] with p the liquid density and R the 

bubble wall velocity, it follows that At = Rc^o/R = 0.9 ps 
(sample A), 0.6 ps (sample B) and ~ 0.3 /is (sample C 
and D). The minimum pressure level which lasts At is 
the negative pressure amplitude Pm of interest. Note that 




FIG. 7: Cavitation bubbles emerging from 6x6 cylindri- 
cal pits with Tc = 246 nm (sample B), for three succes- 
sively applied pressure pulses: a) pm = — 0.24MPa; b) 
= -0.35 MPa; c) pm = -0.54 MPa. The full pattern 
develops when pm = —0.54 MPa is applied without the other 
two preceding pulses (d). 



the difference with the absolute minimum pressure level 
is in most cases only a few percent. 

Let us first consider the case of sample B (rc = 246 nm, 
dc = 2rc). The sample was put 25 mm out of the acoustic 
focus, and three pressure pulses with increasing strength 
(Fig. [4|) were applied successively, without taking the 
sample out of the water. A typical result is shown in 
Fig. [71 With the first pulse {pm = -0.24 MPa) no cavi- 
tation bubbles could be observed (Fig. [7^). The second 
pressure pulse {pm = —0.35 MPa), resulted in a few nu- 
cleated bubbles, but the majority of the nuclei in the 
pattern still did not cavitate (Fig. [Tb). The third pres- 
sure pulse, with Pm = —0.54 MPa, was able to nucleate 
all remaining nuclei (Fig. [7b). Note that the nuclei which 
nucleated with the second pulse could not be nucleated 
with the stronger third pulse, implying that the nuclei 
can be used only once. Nuclei which were not nucleated 
during the second pulse, however, survived and were nu- 
cleated with the stronger third pulse. The experiment 
was repeated eight times with reproducible results (Ta- 
ble [Tj): on average 0, 10, and 26 bubbles were counted 
for the three applied pressure pulses, respectively. When 
the third pulse {pm = —0.54 MPa) was applied without 
the other two preceding ones, the full pattern became 
visible (Fig. [TJl). Since the majority of the pits nucle- 
ated at Pm = —0.54 MPa it is reasonable to assume that 
this pressure amplitude is above the experimental nu- 
cleation threshold for sample B, while the weaker pulse 
{Pm = —0.35 MPa) is (just) below the experimental nu- 
cleation threshold. 
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FIG. 8: Cavitation bubbles emerging from 6x6 cylindri- 
cal pits with rc = 495 nm (sample A), for three succes- 
sively applied pressure pulses: a) pm = —0.20 MPa; b) 
= -0.23 MPa; c) pm = -0.34 MPa. The fuh pattern 
develops when pm = —0.34 MPa is applied without the other 
two preceding pulses (d). 



A similar experiment was carried out with sample A 
{tc = 495 nm, dc = 2rc^ Fig. [8|). Again the sample 
was subjected to three successive pressure signals of de- 
creasing negative pressure without being taken out of 
the water. For the lowest pressure amplitude {pm = 
—0.20 MPa) no cavitation bubbles could be detected op- 
tically (Fig. [8^). A larger amplitude oi Pm = -0.23 MPa 
yielded 14 bubbles of different sizes, with some of them 
barely visible (Fig. [8)3), while a further reduction of the 
liquid pressure {pm = —0.34 MPa) resulted in no visible 
bubbles at all (Fig.[8l3). What happened with the remain- 
ing 36 — 14 = 22 pits? As the lowest negative pressure 
was not able to nucleate them, they must already have 
been nucleated during the first two pulses, i.e. the nucle- 
ation took place below optical resolution. This is possible 
as the resolution of our optical detection is limited and 
the pressure pulse in this case is relatively weak (i.e. R 
is small). Therefore, in contrast with case B, we are not 
able to measure the pressure for which nucleation does 
not take place. When the third pulse {pm = —0.34 MPa) 
was applied without the other two preceding pulses, the 
full pattern became visible (Fig. [8]i). Hence, this is the 
pressure level for which we are sure that full nucleation 
takes place. 

Finally, sample D was studied, which consists of 
nanopits with 50 nm < Tc < 60 nm and varying depths. 
Just like sample C, the shallow pits on sample D nucle- 
ated only a few times, and could not be nucleated in later 
experiments. Fortunately, two columns with the deepest 
pits {dc = 500 nm h Vc = 57 nm, and dc =1000 nm & 
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Amount of bubbles 
3 4 5 


in exp. 1-^ 
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mean 


-0.24 





























-0.35 


11 
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10 
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7 


10 


-0.54 


23 


27 


28 


19 


25 


27 


27 


29 


26 



TABLE I: Results of eight experiments with sample B (rc = 246 nm, dc = 2rc). In each experiment the minimum pressure 
is decreased in three successive steps. The majority of the pits nucleates at pm = —0.54 MPa. 




FIG. 9: Cavitation bubbles emerging from 2x6 cylindri- 
cal pits (sample D) with rc = 57 nm, dc = 500 nm (left col- 
umn) and Tc = 60 nm, dc = 1000 nm (right column), for three 
successively applied pressure pulses: a) pm = —2.3 MPa; b) 
-2.6 MPa; c) pm -3.0 MPa. 



Tc = 60 nm) could be nucleated repeatedly, and the nu- 
cleation threshold could be measured for these pits. The 
experiment was very similar to the ones described before, 
but now the sample was 12 mm away from the acous- 
tic focus. A typical experimental result is depicted in 
Fig. [9l First, a pressure pulse with pm = —2.3 MPa was 
applied and 1 pit from the right column (rc = 60 nm) was 
nucleated. A stronger second pulse {pm = —2.6 MPa) 
was able to nucleate the remaining 5 pits from this col- 
umn, though other pits in the sample did not nucleate, 
as they were smaller. Reducing the negative pressure 
further {pm = —3.0 MPa) resulted in the nucleation of 
the left column of pits with Vc = 57 nm. Hence, a small 
variation in pit sizes of just a few nm is reflected in a dif- 
ferent nucleation threshold. It is also observed that the 
pits did not nucleate a second time, despite their huge 
aspect ratios. 



C. Comparison with theoretical prediction 

How do the experimental results compare with theo- 
retical predictions? In Fig. [10] the theoretical nucleation 
threshold (line), based on Eq. [21 is plotted as a func- 
tion of the pit radius Tc together with the experimental 
results (symbols). We used dc = 2rc as is the case in 
sample A and B. Note that for the pits present in sam- 
ple D the gas pressure term is negligible: the difference 
between dc = 2rc and dc = 20rc changes the theoretical 
prediction for pits of Vc = 50 nm < 1%. The experimen- 
tal data points at which full nucleation was detected for 
samples A, B and D are depicted by crosses, while the 








100 200 300 400 500 
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FIG. 10: Nucleation threshold as function of the pit radius 
for both theory (line) and experiment (symbols, crosses: nu- 
cleation, circles: no nucleation). The inset shows a zoom in 
with errorbars. For visibility overlapping points are shifted 
zbO.25 nm with respect to each other. 

experimental pressures where full nucleation (just) did 
not happen are marked with circles. The inset shows the 
experimental results for sample D including typical er- 
ror bars depicting the standard deviation of the pressure 
recordings. 

We observe a striking quantitative agreement between 
theory and experiment for all samples. Pressure ampli- 
tudes for which nucleation was first detected are below 
the line marking the nucleation threshold. Pressure am- 
plitudes for which nucleation did not occur are either 
above this line, i.e. in the regime where nucleation is not 
expected, or the line is within experimental errorbars. 
For sample D it was observed that the nucleation thresh- 
old strongly depends on the size of the pits: pit radii 
just a few nm smaller resulted in a significantly lower 
nucleation threshold, in agreement with the steep slope 
of Pm(^c) around these values. 



D. Deactivation of cavitation nuclei 

It is well known that artificial nucleation sites in boil- 
ing continue to be active for a long time, emitting many 
bubbles [29]. Similarly, the microscopic wall cracks and 
scratches in a glass full of beer or champagne are seen 
to emit bubbles for a very long time [30^]. Even in cav- 
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it at ion studies on bare hydrophobic substrates, bubbles 
trapped in locahzed defects could be nucleated more than 
a hundred times [19]. Thus, there is something special in 
the deactivation of nuclei observed here which makes this 
situation different from the others. What is the physical 
mechanism responsible for the deactivation of the nuclei? 

The main mechanism responsible for the deactivation 
of nuclei is encountered during the collapse phase of the 
bubble, where a wall-directed jet is formed. In the case of 
a single bubble (or weak interaction among bubbles) the 
jet momentum is expected to be directed mostly normal 
to the wall j3l|, [32[. When bubble-bubble interaction 
is non- negligible, on the other hand, the jet momentum 
may be deflected away from the normal [21]. 

In order to shed light on this proposed deactivation 
mechanism, numerical simulations were carried out to 
elucidate the shape of the air-liquid interface during the 
bubble collapse. For this we used the boundary-integral 
method described in Ref. [33] based on a potential flow 
description of the liquid dynamics. The liquid-solid an- 
gle was prescribed to be Oa = 124° which corresponds 
to the experimentally determined advancing contact an- 
gle of water on the substrates. At the starting point 
of the simulations the bubble was assumed to be a seg- 
ment of a sphere with a radius significantly larger than 
the cavity radius Tc, see Fig. [TT](a). The pressure inside 
the bubble was assumed to be uniform, satisfying the 
adiabatic relation Po/Vq'^ = Pbub/V^^^ with the initial 
pressure po and Vq the volume of the crevice, and p^ub 
and Vbub the instantaneous pressure and volume of the 
bubble, respectively. After release the bubble begins to 
shrink rapidly due to both surface tension and the low 
internal pressure. Eventually it evolves into an almost 
cylindrical shape as illustrated in Fig. [TT] (c). This air 
cylinder collapses radially and finally closes in a single 
point on the axis of symmetry, leaving a small air bub- 
ble entrapped above the pinch-off point, see Fig. [TT] (d). 
Towards pinch-off the liquid rushing radially inwards has 
to accelerate more and more to satisfy the requirement 
of mass conservation. When the advancing liquid front 
reaches the axis of symmetry a high pressure develops 
and the flow is deflected up and down to form two fast, 
needle-like water jets. The continuing collapse of the air 
cavity below and above the pinch-off point provides addi- 
tional momentum to the two jets [34]. The downward jet 
protrudes deeply into the cavity until it hits the bottom 
of the cavity as illustrated in Fig. [TT] (e)-(f). For sim- 
plicity we neglected the upper bubble here which is ex- 
pected to have only negligible influence on the downward 
jet. Upon reaching the bottom the impacting jet would 
form a violent, non-axisymmetric splash which cannot be 
captured by our numerical technique. Nevertheless, one 
can easily imagine the continuation of the process: as 
more and more liquid enters the cavity through the jet, 
the cavity is flooded with liquid, making a second nu- 
cleation impossible. We checked that the jet mechanism 
is present for pits in the size range studied here (50 nm 
< Tc <500nm) and is independent of the initial bubble 
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FIG. 11: (a) Initial configuration for a spherical bubble (blue) 
with radius 5rc on top of the cavity (red). The (advancing) 
contact angle is ^ = 124°. Due to the low pressure inside the 
cavity the bubble starts to collapse (b) evolving into an almost 
cylindrical shape (c), which eventually closes on the axis of 
symmetry in a single point (d). From the pinch-off location a 
downward jet protrudes into the cavity (e) eventually hitting 
the cavity bottom (f). Here it would cause a splash filling the 
cavity with liquid. 



size, pit depth, and contact angle. 

The previous explanation is not applicable to the 
case of strong mutual interaction between the bubbles, 
when the jet tends to be deflected away from the wall- 
normal [21], although the nuclei are still observed to be 
de-activated after emission of the first bubble. However, 
during the expansion phase of the cavitation bubble (typ- 
ically tb ~ 10/is), gas diffuses from the cavity into the 
bubble, thus leaving the pit. The typical diffusion length 
scale is ^ \^Dtb = 10/im, using the diffusion coefficient 
D - IQ-^mVs. In the cases studied here the pits are 
< 1/im deep, i.e. much smaller than the typical diffu- 
sion length scale. Also, the volume of a typical bubble 
is 10^ — 10^ larger than the volume of the nanopits, al- 
lowing the majority of the gas to move from the pit into 
the bubble. [40] Unlike the wall- normal case the different 
jet dynamics does not return the gas to the pit, which 
therefore remains full of vapor and is easily filled by the 
liquid. These conclusions agree with our experimental 
results. In the cases A-D we observed no differences be- 
tween strong (Fig. [6]) and weak (Figs. [7J [8] & [9]) bubble 
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FIG. 12: A super hydrophobic pit (left) can be nucleated hun- 
dreds of times, provided that the liquid jet in the bubble col- 
lapse phase is not directed into the pit. Center: a hexago- 
nal pattern of superhydrophobic pits (100 /xm in between the 
pits) after 230 nucleation events shows only 2 defects. Right: 
a square pattern (300 /am in between the pits) is completely 
intact after 100 shots. 

interaction: in both situations the pits were emptied af- 
ter one nucleation event, in line with the explanations 
provided here. 

This situation can be compared with the previously 
mentioned continuous and long-lived emission of bubbles 
in carbonated beverages and the related phenomena ob- 
served with enhanced surfaces frequently used in boiling 
heat transfer [35|, [s^. In none of these cases the bub- 
bles collapse and therefore the primary mechanism for 
the filling up and consequent de-activation of the pits is 
present. 

E. Superhydrophobic nuclei 

To show the importance of the liquid jet into the pits, 
experiments using superhydrophobic nuclei are illustra- 
tive. Cylindrical pits, etched in Si(lOO) with a diameter 
of 4 /im, were created with a superhydrophobic bottom 
layer (see Fig. [12]), consisting of hydrophobic pillars of 
~ 100 nm in diameter (created through a black silicon 
etching process [37] and a hydrophobic fluoro-carbon 
top layer). The combination of hydrophobicity with 
roughness is known to create superhydrophobicity [sH 
with typical contact angles > 160°. When the bubbles 
were strongly interacting (leading to wall-parallel jets), 
we observed that the pits remained active nucleation 
sites even after hundreds of shots (see Fig. [T2|) . On the 
contrary, in the case of a single bubble (wall-normal 
jet) we observed that the micropits were deactivated 
after a few nucleation events. To explain this striking 
difference one really has to take the direction of the 
liquid jet into account. Apparently, the wall-parallel jet 
is not able to wet the superhydrophobic bottom of the 
pits, while the vigorous, ultra-thin jet directed towards 
the superhydrophobic bottom layer presumably pushes 
the liquid from the dewetted into the wetted state. 



From other work it is indeed known that a force may be 
required to overcome the energy barrier associated with 
this wetting transition [s^. Once in the wetted state, the 
superhydrophobic pit is deactivated and cannot be nucle- 
ated again, apart of course from drying the whole sample. 

V. CONCLUSION 

In conclusion, Atchley & Propseretti's 1989 crevice 
model of cavitation nuclei is experimentally verified using 
nanoscopic well-defined nuclei. Advanced etching tech- 
niques allowed us to create cylindrical pits down to 50 nm 
in radius with high accuracy in both their size, depth 
and mutual position. Upon immersion in water, the hy- 
drophobic nanopits trapped air and served as nucleation 
sites. Stepwise lowering of the acoustic minimum pres- 
sure allowed us to determine the nucleation threshold at 
which the pits start to cavitate. We found that the ex- 
perimental results are in very good agreement with the 
theoretical predictions. This implies that in shock wave 
experiments the size of cavitation nuclei can be deter- 
mined by measuring the pressure at which they start to 
nucleate. 

Cavitation nuclei were deactivated after a single nu- 
cleation event, despite differences in width, depth and 
aspect ratios of the pits. The two mechanisms contribut- 
ing to this effect are diffusion of gas out of the pit during 
the lifetime of the bubble and the subsequent aspherical 
collapse of the bubble. Numerical simulations show that 
in the case of weak bubble-bubble interaction, a sharp 
wall-normal liquid jet is formed which hits the bottom 
of the cavities, thus vigorously wetting the pits. Super- 
hydrophobic nuclei can only be wetted thanks to this 
wall-normal jet. For strong bubble-bubble interactions 
with wall-parallel jets, superhydrophobic pits remain ac- 
tive nucleation sites, even after hundreds of nucleation 
events, in contrast to standard hydrophobic pits. In sys- 
tems were one wants to control the number of cavitation 
nuclei which do not deactivate, superhydrophobic pits 
may find applications. 
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